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We show that two dimensional QCD can, to a good approximation, describe the hadronic structure 

functions measured in Deep Inelastic Scattering. We transform this theory into a new form, Quantum 

HadronDynamics (QHD) , whose semi-classical approximation is closer to nature. The Baryon is 

- - - then a topological soliton, and its structure function can be predicted by a variational principle. This 

G^ 

^\ ' prediction can be tested by comparison with measurements of neutrino scattering cross-sections. 
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Perturbative QCD allows us to determine the dependence of the structure functions of Deep 
Inelastic Scattering ||^. The initial condition for the evolution of the DGLAP equations |2[(i.e., 
, the dependence on the Bjorken scaling variable xb at an initial Q^) cannot be determined within 

00 . the perturbative framework. As a result much effort has been expended to extract these initial 

, distributions by fits to data (See e.g., [^.) We describe in this paper how to derive the valence 

. quark distributions of the proton ( in particular the dependence on a;s at a low initial Q^) from first 

. principles of QCD, by a succession of approximations. It is possible to test the predictions against 



experimental data: a comparison with the measurement of the structure function xFs in neutrino 
scattering will be done in a companion paper The agreement is quite good, which confirms 
the theoretical framework advocated in this paper. There are also some other approaches to this 



^ ' problem which are more numerical in character 

As we will explain below. Deep Inelastic Scattering can be explained by the dimensional reduction 
■ of QCD to two spacetime dimensions. The main idea is now to rewrite two dimensional QCD in terms 

$_( ' of operators that describe mesons rather than quarks and gluons [0. This new formulation which 



I called Quantum HadronDynamics (QHD) has the advantage that its semi-classical approximation 
is quite close to nature: it corresponds to the large Nc limit of QCD. Recall that the semi-classical 
approximation to QCD itself is invalid except at short distances. In particular it fails to explain the 
formation of hadrons. 

In QHD, the baryon appears as a toplogical soliton. Its structure functions are determined 
by a variational principle, within the large Nc limit. There is a natural variational ansatz which 
corresponds to the valence quark approximation. Within this ansatz we can even take care of the 
leading effect of A'^c being finite: it just amounts to restricting the range of momenta allowed for 
partons. Thus we will be able to obtain a variational principle for the valence quark distribution 
functions. A more detailed version of this argument can be found in |8|. 
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Let us begin by recalling why Deep Inelastic Scattering can be understood within a two dimen- 
sional framework [^]. An electron (or neutrino) scatters against a proton (or other nucleus) with a 
space like momentum transfer of magnitude q, causing the target to disintegrate into many hadrons. 
The two Lorentz invariant variables that describe this process can be chosen to be = —q^ and 

2 

the Bjorken variable xb = ~2p^- easy to check that < xb < 1- Deep inelastic scattering is 
the limiting case ^ oo keeping xb fixed. More precisely Q >> where a is the characteristic 
size of the target. In this limit, in the center of mass frame the target will look 'flattened out' 
to a pancake shape due to Lorentz contraction. To first approximation, it can be thought of as 
having infinite extent in the two tranverse directions while of finite extent a in the longitudinal 
direction. In other words the momenta of the constituents of the hadron can be taken to be zero 
in the transverse direction, the corresponding fields are then independent of the transverse spatial 
co-ordinates. Thus the theory of strong interactions (QCD) can be dimensionally reduced to 1 -f 1 
space-time dimensions in describing deep inelastic scattering. We will use the methods of Ref. ^ to 
study this two dimensional field theory. This will give us the structure function as a function of xb- 

The effects of the transverse momenta can then be included as a perturbative correction. As em- 
phasized by Altarelli and Parisi this is precisely the meaning of the DGLAP evolution equations, 
which give the dependence of the structure functions on Q^. (The upper cut-off on the allowed 
transverse momentum is Q, which is how it is related to transverse momentum effects.) Since this 
part of the story is standard [Q] we will not discuss it here. 

In the above two dimensional approximation, the transversely polarized gluons will appear as a 
pair of scalar fields. Since we don't need the gluon structure functions for now, we can ignore these 
transverse gluons. To leading order, the evolution of the valence parton distributions decouple from 
the gluon distribution functions. The longitudinal component of the gluon cannot be ignored as it 
is responsible for the binding of the quarks into the hadron. However, they don't have dynamical 
degrees of freedom and can be eliminated using their equations of motion. 

Thus the action of the two dimensional field theory describing the Deep Inelastic Structure func- 
tion is. 
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S = ^ / tr F^^F'^-d^ + hn ■ V + majqac^d^'x. (1) 

a — 1 

Here, a — 1, • • • , TVc is the color index and a — 1, ■ • • , 2Nf a fiavor index. Also, g is a two dimensional 
Dirac spinor ; each four dimensional spinor splits into a pair of these upon dimensional reduction 
and hence the two dimensional theory has twice the number of flavors as the four dimensional 
theory. In null co-ordinates (for more about the kinematics in the null co-ordinates, see [^) ,and in 
the light-cone gauge A- — 0, we can eliminate the remaining gluon degrees of freedom to get the 
hamiltonian 

J dxx^"^[p+^]xa^ - J ^\x - y\ : X^''\x)xaj{x) :: X^''Hy)XM{y) : dxdy (2) 

where q = I ^ 1 . The upper component qi is not a propagating degree of freedom and has 
\ '72^ ' 

been eliminated in terms of x- The quark fleld x satisfles the canonical anti-commutation relations 
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[x{x),xiy)]+ ~ S{x — y), [x{^) J {y)]+ = 0- Also, the normal ordered product : AB : is defined 
with respect to the vacuum X^(p)|0 for p < 0, x{p)\0 >= for p > 0. 

Now define the color invariant variable M^(x,y) — : [xbaix),x^'^°'{y)] '■ which can be thought 
of as the field operator for a meson field. The space-time points x,y lie along a null line which is 
thought of as the initial value surface. 

Now the entire theory can be described in terms of this color invariant variable. Within the 
subspace of color invariant states, A4!^{x,y) is a complete set of observables; the only operators that 
commute with them are multiples of the identity. This follows from the fact that M{x,y) provide 
an irreducible (projective) unitary representation of the infinite dimensional unitary Lie algebra: 

{ill{p,q), M^r, s)} = ^ I^5l2^5(q - rWd sgn (P - s) + Ma{p, s)] - 5S27r5(s - p)[5l sgn {r - q) + £ll{r, q)]^ . (3) 

Here Mi,{p,q) — J Ml^{x,y)e^''^~"'^dxdy. Note that the commutators are of order so that the 
large Nc limit is a sort of classical limit: plays the role of h in an ordinary field theory. 

In this classical limit the above commutators tend to Poisson brackets of a set of classical dynam- 
ical variables M^{x,y). The phase space of this classical dynamical system must be a homogenous 
symplectic manifold, this being the analogue of an irreducible unitary representation. From the stan- 
dard Kirillov theory (adapted to infinite dimensions by Segal, this phase space is a co-adjoint 
orbit of the unitary group, the Grassmannian. It is the set of all inifinite dimensional operators M 
with integral kernel M^(x,y) satisfying the nonlinear constraint [e-f- A/]^ = 1. Here, e is the Hilbert 
transform operator, €{x,y) = J e'f'^"?') sgn (p)|f . It is also possible to verify the identity above 
directly on color singlet states as is shown in the appendix to Ref. Thus in the large Na limit, 
our problem reduces to solving the equations of motion obtained from the hamiltonian 

E[M] = ~\ [[p+ll]M{p,p)± + 9l [ M^{x,y)Ml{y,x)\x - y\dxdy (4) 
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with the Poisson brackets 

^{M^{x,y), M2{z,u)} = 5l5{y - z)[e'i{x,u) + M2{x,u)] - ^^^(x - u)[el{z,y) + M^,{z,y)]. (5) 

The parameter is related to the quark masses by a finite renormalization: = m\ — where 
nia is the current quark mass. (Also, M^{p, q) — J Ai^{x, y)e~^^^^^''^ dxdy is the Fourier transform.) 

What kind of solution to this theory represents the baryon? The quantity B = J M2{x,x)dx 
can be shown to be an integer, a topological invariant [^. From the definition of M in terms of x, X^ 
we can see that this is in fact the baryon number. Thus the baryon is a topological soliton in this pic- 
ture: an idea originally proposed by Skyrme in quite a different context, and revived by Balachandran 
et. al. and by Witten et. al. We seek a static solution (minimum of the energy subject to con- 
straints) that has baryon number one. Again from the definition in terms of the quark fields, we can 
see that — (p, p) represents the quark number density in momentum space,when p > 0. Similarly, 
^Ma{~p, —p) represents the anti-quark number density. It is convenient to assume a variational 
ansatz of the separable (rank one) form (p, g) = —2il)°'{p)^l}l{q). This satisfies the constraint if 1/; 
is of norm one and of positive momentum: 'ipa{p)\^^ = 1, V'a(p) = 0, for p < 0. This vari- 

able satisfies the Poisson bracket relations {i/)a(p), i^b(g)} — 0, {''/'^(p), ^**(9)} = ~j27r<5*(5(p — q). 
The energy becomes then, 
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a ^ ab 

This variational ansatz corresponds to the valence quark approximation: the anti-quark distributions 
are identicaly zero. In forthcoming papers with V. John and G. S. Krishnaswami, we will show that 
in the limit of zero current quark mass the exact minimum of the energy functional is of this 
separable form In fact deviations are small even for finite m; in the language of the parton 

model, there is a less than one percent probability of finding an anti-quark in the proton at low Q^. 
The factorized ansatz amounts to ignoring the anti-quarks. 

So far we have worked in the large Nc limit. What is the first order effect of A'^c being finite? If we 
stay within the valence parton approximation as we make Nc finite this can be studied by replacing 
the Poisson brackets of i/) by canonical commutation relations: 

As usual, classical Poisson brackets go over to quantum commutation relations, except that the 
role of h is played by These commutation relations have a simple representation in terms of 
bosonic creation-annihilation operators. The constraint on then becomes the condition that there 
be exactly A'^c such bosons in any allowed state. The large Nc limit is then like a thermodynamic 
limit, in the canonical ensemble. 

What are these bosons? A moment's reflection will show that they are in fact the valence partons; 
we have just given a derivation of the valence parton model from a series of approximations on 
QCD. They behave like bosons (rather than fermions) because we are not counting explicitly the 
color quantum rmmber. The wavefunction of the system is completely anti-symmetric in color (to 
make it color invariant) so that the Pauli principle requires it to be symmetric in the remaining 
quantum numbers. The null momentum of each parton is positive so each has to be less than 
the total momentum P. The main effect of a finite (but large) Nc is thus to require that p < P: 
ipijpi, ■ ■ ■ ,Pn) ~ 0, if Pi > P which is in addition to the condition that pi > 0. This way we 
derive exactly the model of interacting partons as an approximation to two-dimensional QCD. 
As noted in that paper, we can make a mean field approximation keeping this condition in place to 
take into account of the effect of finite Nc- 

Thus we can determine the valence parton distribution function if we can minimize the above 
energy functional -Ei(i/') subject to the conditions that 



V(p) = Oforp<Oandp>P, J |^(p)l^^ = 1, J p\^(p)\'^ ^ fp, (g) 

Here / is the fraction of the momentum carried by the valence partons. 

This problem has been solved numerically as well as in a variational approximation . The 
results can then be compared to the experimental measurements of the xF^ structure function. The 
agreement is quite good, confirming our picture of the structure of a hadron In particular we 
have resolved the apparent difference between the Skyrme model of the baryon and the valence 
parton model: we have found a topological soliton model that applies to high energy scattering from 
which an interacting valence parton model can be derived. We can derive more information such as 



4 



spin and flavor-dependent or anti-quark distributions functions and we hope to address these 
issues in later papers. 
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